subroutine fmi_t1ei_fem(fmi, t1, fme, icore, fact)
!
use mod_orbit
use mod_ioff 
implicit none
!
real*8, intent(in) :: fact
real*4, intent(in) :: fmi(*)
real*4, intent(in) :: fme(*)
!
real*4, intent(inout) :: t1(*), icore(*)
!
integer :: i, j, a, b, nsize
real*8, external :: dnrm2
! one should note that fiai in this part actually equal to -fiai
!this is to do fmi = fmi + 0.5*t1(e,i)*fme(e,m)
 call VmnpqVmnrs_to_Vpqrs(fme(ifvo(1)), t1(ifvo(1)), 1, 1, vrta, zorb, 1, & 
                  popa, zorb, 1, popa, zorb, 1, fmi(ifoo(1)), 1, fact, 1.d0)
!
return
end
